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Abstract 

A  short  and  direct  proof  is  given  for  the  convexity  of  the  largest  singiJar  value  of 
e^  Aoe~  on  any  convex  set  of  commuting  matrices  in  i2"^",  where  Ao  €  R^^^  is 
fixed.  This  statement  does  not  hold  if  X  is  permitted  to  be  a  general  square  matrix. 
A  counterexample  is  supplied  for  noncommuting  matrices. 


The  work  of  both  authors  was  supported  in  part  by  National  Science  Foundation 
Grant  DCR-85-02014. 


Doyle's  structured  singular  value  analysis  [l]  is  a  well  known  tool  in  robust  control. 
One  technique  (see  [1,2])  used  to  bound  the  structured  singular  value  is  to  minimize 
the  largest  singular  value,  i.e.  the  two-norna,  of  the  matrix  DAoD~^ ,  where  Aq  is 
fixed  and  D  is  a  variable  positive  diagonal  matrix.  Doyle  and  Chu  showed  that 
the  largest  singular  value  of  this  matrix  is  a  convex  function  of  X,  where  D  =  e  . 
Hence,  the  optimization  problem  is  tractable.  Their  proof  consists  of  two  steps, 
neither  of  which  is  readily  accessible  in  the  open  literature.  The  first  step,  a  char- 
acterization of  convexity,  may  be  found  in  [3]  and  the  second  step  may  be  found  in 
[4]. 

In  this  note  we  give  a  new,  short,  and  direct  proof  of  the  convexity  property.  We 
show  that  the  theorem  applies  to  any  convex,  commuting  set  of  matrices  X .  We 
also  show  that  the  result  does  not  hold  if  X  is  permitted  to  be  a  general  square 
matrix. 

Theorem:     Let  fl  be  a  convex  subset  of  i?"''",  such  that  X,Y  e  ^  =J>  XY  =  YX. 
For  any  Aq  G  iZ"'*",  the  l2-noTux  of  the  matrix  valued  function 
A-.n^  iZ"''"  defined  by 

A{X)  =  e^  Ao  e-^  (l) 

is  convex  on  CI. 

Proof:     Since  ||>i(X)||  is  a  continuous  function,  it  suffices  to  show  that 

-\\A{X)\\  +  -\\A{Y)\\  >  \\A{Z)\\  ,      where  Z  =  {X  +  Y)/2,   WX,Y  e  Q,      (2) 

2  2 

in  order  to  prove  convexity  [5,  p.  101].  In  the  following  argument,  5""^  denotes 
the  unit  sphere  {^  e  R^\    \\i\\  =  1}-  We  define  the  matrix 

£;^,(y-X)/2^  (3) 

and  the  right  and  left  singular  vectors,  v  and  u,  of  A{Z)  corresponding  to  its  largest 
singular  value: 

A{Z)v  =  \\A{Z)\\u  (4) 

A^{Z)u  =  \\A{Z)\\v  (5) 


Using  the  triangle  inequedity  and  the  hypothesis  XY  =  YX,  we  have: 


-\\A{X)\\  + -\\A{Y)\\       > 

2  2 

>  '-\\A{X)  +  A{Y)\\ 

=  -\\E-'^A{Z)E  +  EA{Z)E 


-1 


sup  C'^{E-^A{Z)E  +  EA{Z)E-'^)( 


>  1      '''^^      {E-'A{Z)E  +  EA{Z)E-')        ^"^ 


2   VuTJEhi 

1     v'^E^v  +  u'^E'^u 


Vv^E^v 


2    Vu^E^u  Vv'^E^v 


\\A{Z)\\ 


v'^E^v  u^E^u 

+ 


uTE^u 


v'^E^v 


\\A{Z)\\ 


>  \\A{Z)\\ 


The  last  inequality  follows  from  the  fact  that 


\H)>-^' 


Vd>0 


(6) 
(7) 
(8) 

(9) 
(10) 

(11) 
(12) 


(13) 


An  example  of  a  convex,  commuting  subset  of  iZ"**",  is: 

*M  =  {m  diag(x)M-^  I     I  G  iZ"} 


(14) 


where  M  G  i?"^",  det(M)  j^  0,  and  diag(x)  denotes  the  n  x  n  diagonal  matrix 
which  satisfies  (diag(i))^^  =  Xi.  In  particular,  ^/  is  the  set  of  diagonal  matrices. 

It  has  generally  been  thought  that  Doyle  and  Chu's  proof  applies  to  general  square 
matrices.  However,  this  is  not  true,  as  the  following  counterexEimple  shows: 


Ao  = 


X  = 


-0.88 

1.00  ■ 

1.00 

0.88 

1.53 

0.22" 

0.22 

4.40 

y  = 

'1.28 
2.31 

2.31  ■ 
5.32 

11^(^)11  +  M^)II~  9.1361 
2                  2 

V  2  ; 

~  20.8545 

(15) 


Finally,  we  note  that  although  ||A(X)||  is  convex  on  fi,  it  is  not  differentiable  ev- 
erywhere, since  the  singvdar  values  of  a  matrix  are  generally  not  differentiable  at 
points  where  they  cotilesce  to  multiple  values.  This  is  well  known  and  treated  in 
[l];  another  relevant  reference  is  [6].  See  [7]  for  a  recent  contribution  showing  how 
to  obtain  an  algorithm  which  is  quadratically  convergent  even  to  nondifferentiable 
solutions.  The  aJgorithm  is  developed  for  the  related  problem  of  minimizing  the 
maximum  eigenvalue  of  a  linear  combination  of  symmetric  matrices;  we  are  cur- 
rently extending  this  to  minimize  (1). 
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